Optical phase-space reconstruction of mirror motion at the attometer level 
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We describe an experiment in which the quadratures of the position of an harmonically-bound 
mirror are observed at the attometer level. We have studied the Brownian motion of the mirror, 
both in the free regime and in the cold-damped regime when an external viscous force is applied 
by radiation pressure. We have also studied the thermal-noise squeezing when the external force is 
parametrically modulated. We have observed both the 50 % theoretical limit of squeezing at low 
gain and the parametric oscillation of the mirror for a large gain. 
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I. INTRODUCTION 

Optomechanical coupling, that is the cross-coupling 
between the motion of a mirror and a laser field reflected 
upon it, first appeared in the context of interferometric 
gravitational-wave detection 0, with the existence of 
the so-called Standard Quantum Limit H EL 01 The 
unique sensitivity of interferometric techniques has since 
then been used for other high-sensitivity measurements, 
such as AFM [|| or optical transducers 0,@. The field 
has recently encountered much attention and has gained 
a life of its own in the quantum optics community: sev- 
eral schemes involving a cavity with a movable mirror 
have been proposed either to create non-classical states 
of both the radiation field 0, 0] and of the motion of 
the mirror pHIT^.[l^ |. or to perform QND measurements 
[Tij . Recent progress in low-noise laser sources and low- 
loss mirrors have made the field experimentally accessible 
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A recent trend in quantum optics is to fully reconstruct 
the quantum state of either a mode of the radiation field 
through the quantum tomography technique [Tsl Il9l | or a 
trapped atom [20|. but to our knowledge, no such exper- 
iment has yet been performed at the quantum level on a 
mechanical oscillator. 

In this paper, we present an experiment which com- 
pletes the analysis of the motion of a plano-convex mir- 
ror given in 0, llH. l2l| , reconstructing the phase-space 
distribution of the motion through the simultaneous clas- 
sical measurement of both quadratures of the mirror po- 
sition. The technique is applied to a variety of states 
of motion: Brownian motion and its cold-damped and 
squeezed counterparts. 

In Sec. [H]we present the experimental setup used to 
monitor the motion of the mirror and we describe how 
the motion in phase-space can be reconstructed. In Sec. 
IIIII we consider the case of a free mirror at thermal equi- 
librium and we compare our experimental results to pre- 



dictions of the fluctuations-dissipation theorem [22j . The 
experimental setup has also been modified to include an 
external viscous force applied to the mirror, and a similar 
analysis is given in Sec. IIVI in the corresponding case of 
the cold-damped regime. Our results demonstrate that a 
new equilibrium is obtained with a reduction of thermal 
noise. 

In Sec. E|we present an experiment of parametric am- 
plification of the thermal noise |2^| , below the oscillation 
threshold. The experimental setup is modified in order to 
modulate the strength of the viscous force, and squeezing 
of the thermal noise has been observed. The observation 
of parametric oscillations of the mirror is finally described 
in Sec. EH 
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II. EVOLUTION IN PHASE-SPACE 

The mirror motion is monitored by an optomechanical 
displacement sensor. It relies on the sensitivity of the 
phase of a reflected light beam to mirror displacements. 
Monitoring such a phase-shift allows to reconstruct the 
dynamical evolution of the mirror. This phase-shift can 
be induced by various kinds of mirror motions, cither 
external 0, or internal [llj HjJ |2(j . The former is 
important for suspended mirrors since the excitation of 
pendulum modes of the suspension system leads to global 
displacements of the mirror. The latter corresponds to de- 
formations of the mirror surface due to the excitation of 
internal acoustic modes of the substrate. These various 
degrees of freedom have however different resonance fre- 
quencies and one can select the mechanical response of a 
particular mode by using a bandpass filter. 

In our experiment we detect only frequencies around 
the fundamental acoustic resonance of the mirror so that 
the mechanical response is mainly ruled by the behaviour 
of this particular internal mode. The mirror motion can 
then be approximated as the one of a single harmonic 
oscillator characterized by its resonance frequency Qmj 
its quality factor Q, and its mass M . This mass actually 
corresponds to an effective mass which describes the ef- 
fective motion of the mirror as seen by the light, that is 
the deformation of the mirror surface averaged over the 
beam spot [MIMI^. 
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The temporal evolution of the mirror position x (i) is 
given by linear response theory [2^|. Assuming that a 
force F is applied to the mirror, the displacement x [Vt\ 
in Fourier space at frequency f2 is related to the force by 



(1) 



where \ [^] is the mechanical susceptibility of the mir- 
ror. For a viscously-damped harmonic oscillator it has a 
Lorentzian behaviour given by 



m (QL - n 2 - im) : 



(2) 



where T = Qm/Q is the damping rate |29| . 

Depending on the nature of the applied force F, the 
spectrum of the mirror motion will typically be peaked 
around the mechanical resonance frequency Qm over a 
frequency span T, which is much smaller than Hm for a 
high-Q harmonic oscillator. To study the temporal evo- 
lution of the motion in phase-space, it is preferable to 
describe the motion in the rotating frame in order to re- 
move the intrinsic oscillatory dependence with time. This 
leads to the two quadratures X\ and X 2 , defined by 

x (t) = Xi (t) cos (fi M i) + X 2 (t) sin (Q M t) . (3) 

This equation yields the expression of the two quadra- 
3r space, 



tures in Fourier space, 



X\ [u>] = x [Qm + w] + x [ 
X 2 [uj] = —i (x [SIm + 10] - 



x [-Q M . 



(4) 
(5) 



Quadratures X\ (t) and X 2 (t) vary very slowly with time 
(over a 1/r timescale), and lo which represents the fre- 
quency mismatch — Qm between the analysis frequency 
and the mechanical resonance frequency is always consid- 
ered small as compared to fljf. 

Expressions of the quadratures in presence of an ap- 
plied force can be derived from eq. Q by inserting the 
expressions of x [£Im + lo] an( i x [— f2j\/ + to] in the defini- 
tions of Xi [lo] and X 2 [lo]. We find after simplification, 
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X 2 M = 
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2MCl M \-ilo + T/2 
1 ( 1 



2MQ M \-ico + T/2 



F 2 [lo] , (6) 
^1 M , (7) 



where we have introduced the two quadratures F± and F 2 
of the applied force, defined in Fourier space by similar 
expressions as for X\ and X 2 (eqs. |Uand|SJ. Both posi- 
tion quadratures have a Lorentzian response centered at 
zero frequency and of width T. 

The harmonic oscillator we consider throughout the 
paper is the high-Q fundamental acoustic mode of a mir- 
ror coated on a plano-convex resonator made of fused 
silica. The resonator is 1.5 mm thick at the center with 
a diameter of 14 mm and a curvature radius of the convex 



side of 100 mm. Such dimensions lead to a resonance fre- 
quency of the fundamental mode of the order of 2 MHz. 
The oscillator's parameters have the following values [l5j . 

tt M = 2vr x 1859 kHz, M = 230 mg, Q = 44000. (8) 

The mirror coated on the flat side of the resonator is used 
as the end mirror of a single-ended Fabry-Perot cavity 
with a Newport high-finesse SuperMirror as input mirror 
(figure The whole provides a 1 mm- long cavity with 
an optical finesse T — 37000. The light beam entering 
the cavity is provided by a frequency-stabilized titane- 
sapphire laser working at A = 810 nm. The light beam is 
also intensity-stabilized and spatially filtered by a mode 
cleaner. 
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FIG. 1: Experimental setup. A plano-convex mirror is used 
as the end mirror of a high-finesse cavity. A frequency and 
intensity-stabilized laser beam is sent into the cavity and the 
phase of the reflected beam is measured by a homodyne de- 
tection. A demodulation system working at the mechanical 
resonance frequency Qm extracts both quadratures of the sig- 
nal. 

Near an optical resonance of the cavity the intracav- 
ity intensity shows an Airy peak when the cavity length 
is scanned through the resonance, and the phase of the 
reflected field is shifted by 27r. The slope of this phase- 
shift strongly depends on the cavity finesse and for a 
displacement 5x of the end mirror one gets at resonance 
a phase-shift dtp on the order of 



Sx 

Sip ~ SJ 7 — + S(p n , 
A 



(9) 



where A is the optical wavelength and 5tp n the phase 
noise of the reflected beam. At high frequency, all tech- 
nical noises can be suppressed and the phase noise corre- 
sponds to the quantum noise of the incident beam which 
is inversely proportional to the square root of the light 
power. 

The phase of the reflected field is monitored by a ho- 
modyne detection working at the quantum level. For a 
100 /j,W incident laser beam, we have shown that the sen- 
sitivity to mirror displacements is limited by the quantum 
noise of light and it has been measured to 



Sx min = 2.8 10~ 19 m/VWz, 
at an analysis frequency of 2 MHz [Tlj . 



(10) 



3 



To extract both quadratures of the mirror motion we 
use a demodulation system at frequency Am (figure |2J). 
A bandpass filter of width 10 kHz first selects frequencies 
around the fundamental resonance frequency in the sig- 
nal given by the homodyne detection. The signal is then 
mixed with two square signals at frequency Qm> ° ne de- 
phased by 7r/2 with respect to the other one. Harmonics 
of the output signals are cancelled out by two electronic 
low-pass filters. The filters also have to cancel the signal 
due to other acoustic modes of both the resonator and the 
coupling mirror. The expected width of the fundamental 
mode being of the order of 43 Hz and the closest mode 
3 kHz above, we use second-order filters with a 500 Hz 
cut-off frequency. We have checked the transfer function 
of each filter in order to ensure that the dissymmetry 
between both channels is less than 1%. 
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FIG. 2: Demodulation of the position signal: the signal tp(t) 
given by the homodyne detection is filtered and mixed with 
two 7r/2-dephased square signals at the mechanical resonance 
frequency Qm. Low-pass filters are used to extract the low- 
frequency component of the quadratures. 

Note that X\ and X2 are conjugate quantum observ- 
ables and therefore cannot be measured simultaneously 
at the quantum level. This is not in contradiction with 
our experimental setup since we extract the position 
quadratures X± and X2 from the phase ip of the reflected 
field which also contains the quantum phase noise 5ip n of 
light (eq. |5J). A measurement of the quadratures is thus 
contaminated by the quantum noise of light so that it 
is not possible to measure simultaneously both quadra- 
tures at the quantum level. For this paper this does not 
constitute a limitation since the signals studied in next 
sections (thermal noise and external force) are way above 
the Standard Quantum LimiMmp_osed by the measure- 
ment and back-action noises 



III. OBSERVATION OF BROWNIAN MOTION 

We present in this section the results obtained when 
the mirror is free (no external force) and at room temper- 
ature. The thermal equilibrium can be described as the 
result of a coupling with a thermal bath via a Langevin 
force Ft applied on the mirror. The resulting Brownian 
motion is the response of the acoustic mode to this force 
according to eq. Q. The fluctuations-dissipation theo- 
rem relates the spectrum St of the Langevin force to the 
dissipative part of the mechanical susceptibility [2^ , 



where ks is the Boltzmann constant and T the temper- 
ature of the thermal bath. 

The sensitivity of our experiment is high enough to de- 
tect the Brownian motion. We have already observed the 
noise spectrum of the mirror displacement, which has a 
Lorentzian shape centered at the mechanical resonance 
frequency VLm and a width T [lfij . The height of the peak 
is at least 4 orders of magnitude larger than the sensitiv- 
ity defined by the quantum phase noise of the reflected 
beam (eq. 1101) . 

We have observed the temporal evolution of the Brow- 
nian motion by sending the outputs of the demodulation 
system in a digital oscilloscope. The resulting traces of 
the two quadratures X\ and X2 are plotted in the upper 
curve of figure for an acquisition time of 500 ms. The 
temporal trajectory in phase-space appears as a random 
walk around the center (X\ = 0, X2 = 0). 




2k B T 

n 



Im(l/xM) 



(11) 



FIG. 3: Brownian motion of the plano-convex mirror in the 
(Xi,X2) phase-space. Full scale corresponds to ±2 10 -15 m. 
Top: temporal acquisition over 500 ms. Bottom: histogram 
of the distribution acquired over 10 minutes. The number of 
cells is 256 x 256, the limit of sensitivity (quantum phase noise 
of light) corresponding approximately to 1 cell. 

This plot can be compared to the Brownian motion 
of a mirror coated on a torsion oscillator, as observed in 
ref. p"?! by a similar optical interferometry technique. 
The larger density of points and the more circular shape 
of our results is due to the larger ratio between the ac- 
quisition time T acq and the evolution time (~ 2-7r/F) of 
the Brownian motion in phase-space. In our experiment, 
Tacq = 500 ms and r/27r = 43 Hz, so that the phase 
space is mapped approximately 20 times. In the case 
of the torsion oscillator, as the mechanical resonance fre- 
quency is two orders of magnitude lower, r/27r ~ 20 mHz 
and the phase space is mapped just once even for a 
Tacq = 80 s acquisition time. This explains why results 
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in ref . |17| look similar to the Brownian motion of a free 
particle: ergodicity is irrelevant on such short timescales. 

Our results are more likely to be compared to the ther- 
momechanical noise observed with an AFM cantilever 
|23| , although the simpler two- wave technique used in this 
work yielded a lower sensitivity, on the order of 10 -12 m. 
In this experiment the resonance frequency also lay in 
the 20 kHz range but the lower mechanical quality fac- 
tor of the cantilever drastically lowered the timescale of 
the Brownian motion, enabling to display results compa- 
rable to ours. 

We have calibrated the observed displacements by us- 
ing a frequency modulation of the laser beam. As shown 
by eq. Q we measure a displacement with a reference 
corresponding to the optical wavelength. A displacement 
Ax of the mirror is thus equivalent to a frequency varia- 
tion Av of the laser related to Ax by 



Av 
v 



Ax 



(12) 



where v is the optical frequency and L the cavity length. 
Plots of figure correspond to a full scale of the os- 
cilloscope of ±100 mV . We have applied a frequency 
modulation of calibrated amplitude Av = 200 Hz and 
obtained a deviation in phase-space of 27 mV. From 
cq. <|12l) this corresponds to an equivalent displacement 
Ax = 5.4 10~ 16 m, so that 1 mV is equivalent to 2 10~ 17 
m and full scale in figure corresponds to ±2 10" 15 m. 

This calibration also allows to determine the sensitivity 
of the measurement in phase-space. The quantum phase 
noise of light induces a rms voltage noise at the output 
of the demodulation system of 0.86 mV. The smallest 
observable displacement in phase-space is then 



Ax„ 



1.7 10~ 17 



(13) 



and corresponds to a few attometers. One can relate this 
sensitivity to the one expressed in term of spectral am- 
plitude (eq. ITU|) . Since the phase noise of the reflected 
beam is a white noise on the frequency scale of the mea- 
surement, the minimum measurable displacement Aa; m i n 
in phase-space is equal to the spectral sensitivity 5x m i n 
integrated over the measurement bandwidth, that is 



Aa4 in = (Sx min ) 2 / \H[w]\ 2 <b* 



(14) 



where H [uj] is the transfer function of the low-pass fil- 
ters at the output of the demodulation system. From the 
value of the spectral sensitivity (eq. HOfl and the char- 
acteristics of the filter (second-order filter with cut-off 
frequency 460 Hz and quality factor 2200), one gets 



Ax„ 



1.65 10 



-17 



(15) 



in excellent agreement with the measured value feci. 1 1 3|> - 
As shown by eq. I|14|l the sensitivity depends on the 
frequency cut-off of the filters and can be increased with 
a smaller measurement bandwidth. 



Acquiring data during a longer time allows to recon- 
struct the distribution in phase-space. The full voltage 
scale at the output of the demodulation system is divided 
into 256 cells so that one cell approximately corresponds 
to the limit of sensitivity Ax m i n . We then perform an 
histogram by accumulating the temporal traces of the 
quadratures in the 256 x 256 cells. For an acquisition 
time of 10 minutes the total number of points delivered 
by the digital oscilloscope is of the order of 6 10 5 and 
the phase-space is mapped approximately 25000 times. 
One then gets a good statistics for the phase-space dis- 
tribution as shown in the bottom curve of figure The 
distribution has a gaussian shape with a revolution sym- 
metry verified with an agreement better than 1%. The 
width of the distribution is 



AX i = AX 2 = 36.6 10 



-17 



(16) 



This result can be compared to the theoretical value 
expected for the Brownian motion of a harmonic oscilla- 
tor. The two quadratures Ft x and Ft 2 of the Langevin 
force are uncorrelated and have a flat spectrum (S 1 ^ = 
St 2 = 2St). From eqs. ©, and (JTTJ one gets the 
spectrum for the two position quadratures 



s Xl M = Sx 2 [uj] 



Mn 2 M (u 2 + r 2 /4) 



k B T. 



(17) 



Both noise spectra are centered around zero frequency 
with a Lorentzian shape of width T. Integration over 
frequency leads to the variances 



AX 2 = AXl = 



k B T 
MQ 2 M 



(18) 



According to the characteristics of the harmonic oscilla- 
tor (eq. |SJ) one gets dispersions for the two quadratures 
of 36.3 10 -17 m in excellent agreement with the experi- 
mental value (eq. I16|) . 

Finally we have determined the correlation function of 
the motion in phase-space defined for quadratures Xi, Xj 
= 1,2) as 



C ij (r)=(X i (t)X J (t + T)) i 



(19) 



where the brackets 



stand for the temporal average 



over the measurement time. Curves a in figure0]show the 
experimental result obtained for an acquisition time of 10 
minutes. As expected for a harmonic oscillator in ther- 
mal equilibrium, there is no cross-correlation between the 
two quadratures (C12 = C%\ = 0) and auto-correlation 
functions C\\ and C22 are equal and exhibit an exponen- 
tial decay from their initial values corresponding to the 
variances, 



C ti (r) = AX 2 exp (-IV/2) 



(20) 



The slope of the curve leads to a damping constant r/27r 
of 48 Hz in good agreement with the expected value 
(r/27r~ 43 Hz). 
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FIG. 4: Plots in lin-log scale of the temporal correlation func- 
tions Cn (r), C22 (t) for both quadratures as a function of the 
delay r. a: free motion at room temperature. 6: cold damped 
mirror (section llVi . c and d: parametric cooling (section FVt. 
e: phase noise of light. 



Curve e in figure 01 represents the correlation function 
for the phase noise of the reflected beam obtained when 
the cavity is out of resonance with the incident laser. One 
gets an exponential decay related to the cut-off frequency 
of the low-pass filters of the demodulation system. This 
curve clearly shows that the measurement noise is almost 
negligible as compared to the thermal noise of the mirror. 



IV. COLD-DAMPED REGIME 

We now study the temporal evolution of the mirror in 
the cold-damped regime obtained by freezing the motion 
with an additional radiation pressure applied on the mir- 
ror. The principle of the cold damping mechanism 30] is 
to monitor the thermal motion and to use an electronic 
feedback loop which corrects the displacement by apply- 
ing an appropriate damping force [16j. 

As shown in figure an acousto-optic modulator and 
an intensity-modulated beam is added to the experimen- 
tal setup. This allows to apply a controlled radiation 
pressure on the mirror. An electronic differentiator with 
a variable gain drives the modulator to create a viscous 
force, that is a force proportional to the measured veloc- 
ity of the mirror. The incident beam on the modulator 
is intensity-stabilized in order to apply a well calibrated 
force. 

The effect of this force is to modify the damping of the 
mirror. It can be described as a change in the mechanical 
susceptibility which now becomes |l6| 



1 



Xfb = 



m (n 2 M - n 2 + (1 + g ) rn) ' 



(21) 



where g is a dimensionless parameter characterizing the 
gain of the feedback loop and proportional to the elec- 
tronic gain and to the intensity of the auxiliary beam. As 
no additional thermal noise is associated with the cooling 
process, the mirror motion is still described by the evo- 
lution equation Q where the mechanical susceptibility 
X is replaced by Xfb- The resulting motion is equivalent 
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FIG. 5: Cold-damping experimental setup. The mirror is 
submitted to the radiation pressure of an auxiliary laser beam 
reflected from the back on the mirror and intensity-modulated 
by an acousto-optic modulator (AOM). The feedback loop 
allows to apply an additional viscous force. 



to a thermal equilibrium but at a different temperature, 
depending on the gain g of the loop. The fluctuations- 
dissipation theorem allows to define an effective temper- 
ature T e ff which is smaller than the room temperature 
T for a positive gain g, 



T eff ~ 



T 



9 



(22) 



The freezing of the mirror of course goes with a reduc- 
tion of the Brownian motion which would appear as a 
shrinking of the distribution in phase-space. The vari- 
ances of quadratures X\ and Xi can be calculated in the 
same way as in previous sections and one gets expressions 
similar to eq. 118J) with temperature T replaced by T e ff, 



AXf = AX? 



M& M (1 



9) 



(23) 



The change in the mechanical susceptibility also affects 
the spectra of X\ and X-i- They still have a Lorentzian 
shape but with a width equal to the effective damping 
(1 + g) r and a power at resonance reduced by a factor 

(1+9) 2 - 

The cold damping of a mirror has already been demon- 
strated by studying the thermal noise spectrum of the 
mirror |2lJ. Noise reductions larger than 30 dB at 
the mechanical resonance frequency and temperature re- 
duction factors larger than 10 have been obtained. We 
present here the results obtained in phase-space for a 
moderate gain g ~ 3. 

The motion of the cold-damped mirror is shown in 
figure HO (top curve) which also displays the related his- 
togram (bottom curve). As expected this distribution 
has a gaussian shape with cylindrical symmetry, but with 
a reduced width as compared to the initial thermal dis- 
tribution of figureOl Ratios of the variance of the quadra- 
tures between the two situations are in agreement with 
cq. <|23[) and with the value of the feedback gain. 

Curves b in figure^Jshow the auto-correlation functions 
for the two quadratures. They correspond to the ones of 
a harmonic oscillator in thermal equilibrium at a lower 
temperature T e ff. Compared to the Brownian motion at 
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room temperature, the initial value at r = is decreased 
as the variances and the time constant of the exponential 
decay is reduced. This is associated with the increase of 
the effective damping by the cold damping mechanism 
((1 + g)T/2-K ~ 170 Hz). 
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FIG. 6: Cold-damped motion with a gain g = 3, in the 
(X\,X2) phase-space (same scale as in figure |3J. The shrink- 
ing of the motion is due to the cooling mechanism. Top: 
temporal acquisition over 500 ms. Bottom: corresponding 
histogram acquired over 10 minutes. 



V. PARAMETRIC COOLING 

Distributions obtained in previous sections are sym- 
metric since there is no privileged quadrature. We now 
explore the possibility to obtain asymmetric distribu- 
tions, that is to squeeze the thermal noise. By anal- 
ogy with quantum squeezed states of light, one conve- 
nient way to achieve such states is to use parametric am- 
plification. One can take advantage of the presence of 
an intensity-controlled auxiliary laser beam to modulate 
some parameters of the mechanical oscillator at twice the 
resonance frequency Q m ■ 

Parametric amplification of a mechanical harmonic os- 
cillator is usually done by modulating its spring constant 
[23L l3l| . An equivalent mechanism is obtained for a mir- 
ror by applying a modulated force proportional to the 
position x of the mirror, 



F (t) = 2gMTQ M cos (2Q M i) x (t) , 



(24) 



where g is the gain of the parametric amplification. This 
can easily be done in our experiment by monitoring the 
mirror position and accordingly controlling the force ap- 
plied by the auxiliary beam. Since we use an electronic 



control it is even possible to differentiate the signal given 
by the homodyne detection as in the case of cold damp- 
ing, so that we apply a modulated force proportional to 
the velocity x of the mirror, 



F (t) = 2gMT cos [2Vt M t) x (t) . 



(25) 



This corresponds to the less usual case of a parametric 
amplification via a modulation of the relaxation rate of 
the mechanical oscillator. 

For such a viscous force one gets from eqs. 0} and (J5J 
the two quadratures F\ and F2 of the applied force, 



F l [u] = gMTQ M X 2 [u] , 
F 2 [uj] = gAimrnXx [uj] . 



(26) 
(27) 



According to the evolution equations (JSJ and J7J one finds 
that the two quadratures X\ and X 2 of the mirror posi- 
tion are decoupled and obey the following equations, 



Xi[uj] = - 
X 2 [lu] 



1 



2M9, M 
1 



1 



Ti/2 



2MQ 



M 



r 2 /2 



Fn 



F Tl M 



(28) 
(29) 



where Ft, and Ft 2 are the two quadratures of the 
Langevin force Ft and where T%, T 2 are the quadrature- 
dependent effective dampings in presence of the modu- 
lated feedback force, 



r x = T(l + g). 
T 2 = T(l-g). 



(30) 
(31) 



As for the cold damping mechanism, the parametric 
amplification changes the damping of the oscillator but 
the effect now depends on the quadrature, one damping 
being increased while the other one is decreased. As a 
result, one quadrature of the motion is amplified whereas 
the other one is attenuated. From eqs. (|28|) and (|29|l one 
gets the variances, 



AAj 2 



AXi 



k B T r 
MSl 2 M Ti ■ 
k B T r 

msi 2 m r 2 ■ 



(32) 
(33) 



which have to be compared to eq. (|18fl . 

In the case of a modulated restoring force (eq. [33J| 
the two quadratures F\ and F 2 of the force respectively 
depend on X\ and X 2 so that the evolution equation of 
X\ and X 2 are no longer decoupled. One has to consider 
new quadratures obtained by a rotation of 45° in phase- 
space in order to obtain decoupled equations similar to 
eqs. (|28H and l|29|). with the same effective dampings Ti 
and T 2 . Both forces thus give the same results except for 
a global rotation in phase-space. 

The feedback loop in our experimental setup is modi- 
fied in order to modulate the amplitude of the force ap- 
plied by the auxiliary laser beam at twice the resonance 
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frequency Am (figure [JJ . A reference signal is used to 
synchronize this modulation with the demodulation of 
the quadratures. We can change the phase between the 
two signals at Q,m and 20, m in order to select the de- 
coupled quadratures, either X\ and X2 for a modulated 
viscous force (eq. 125(1 or quadratures rotated by 45° for 
a modulated restoring force (eq. 
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FIG. 7: Experimental setup of parametric amplification. An 
electronic mixer modulates the signal of the feedback loop at 
twice the mechanical resonance frequency Qm- The demod- 
ulation system is synchronized with the reference signal at 
20m • 




FIG. 8: Parametric amplification of the Brownian motion for 
a gain g = 0.8, in the (X\,X2) phase-space (same scale as 
in figure [3J . The squeezing of the thermal noise is due to 
the quadrature-dependent cooling. Top: temporal acquisition 
over 500 ms. Bottom: corresponding histogram acquired over 
10 minutes. 

The resulting motion in phase space is shown in figure|Hl 
for a parametric gain g — 0.8. The distribution still has 
a gaussian shape but with different widths for the two 
quadratures. The distribution is no longer symmetric 
and dispersions of both quadratures are found equal to 

(34) 



Compared to the dispersion obtained at room tempera- 
ture (eq. [HI the variance of quadrature X\ is reduced 
by a factor 1/ (1 + g) whereas the one of quadrature X 2 
is increased by a factor 1/ (1 — g). One quadrature is 
thus cooled while the other one is heated and the state 
corresponds to a thermal squeezed state. 

Curves c and d of figure 01 show the auto-correlation 
functions for the two quadratures X2 and X\, respec- 
tively. Initial values (r = 0) correspond to the variances 
deduced from eq. (|34|1 and the time constants of the ex- 
ponential decays are related to the effective dampings T2 
for curve c and Ti for curve d, in agreement with eqs. i|30|) 
and l|31|) . Note that we have aligned the two quadratures 
X\ and X2 with the axes of the squeezed thermal state 
by a proper choice of the phase between the two refer- 
ence signals at SIm and 2Q,m- We have then checked that 
there is no cross-correlation between the two quadratures 
(Ci 2 = C21 = 0). 

The treatment presented here is based on linear re- 
sponse theory and is only valid for a gain g smaller than 
1. For a gain equal to 1, the effective damping T2 van- 
ishes fea. !31|) and the variance of quadrature X2 becomes 
infinite (cq. . This actually corresponds to the oscil- 
lation threshold of the parametric amplification and will 
be studied in next section. 




AXx = 27 1CT 17 m, AX 2 = 78 10 
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FIG. 9: Relative values of effective dampings Ti (a), F2 (6), 
and variances AXl ( c) , A-Xf ( d) normalized to their values in 
thermal equilibrium at room temperature, for different gains 
of amplification. Solid lines are theoretical predictions. 

We have checked the theoretical behaviour of the para- 
metric amplification as a function of the gain g in the 
linear regime (g < 1). We have repeated the experiment 
for different values of the gain and reported the effective 
dampings Tj and the variances AXf for the two quadra- 
tures (i = 1,2). Points in figure^lshow the experimental 
results for these parameters, normalized to their respec- 
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tive values in thermal equilibrium at room temperature. 
The gain g is estimated by an average over the indepen- 
dent values of the effective dampings and the widths of 
the distributions in phase-space. 

The measurements are in very good agreement with 
the theoretical expressions given by eqs. I|30[) to (|33[1 and 
shown as solid lines in figure without any adjustable 
parameter. Both dampings have linear and opposite de- 
pendence with the gain. The variance of quadrature X2 
increases and diverges at the oscillation threshold g = 1, 
whereas the one of quadrature X x goes down to 0.5. The 
efficiency of parametric amplification which is mainly re- 
lated to the power of the auxiliary laser beam is large 
enough to obtain a gain g close to unity. One then reaches 
the 50% theoretical limit of squeezing for the variance of 
quadrature X\. 



VI. PARAMETRIC OSCILLATION 

As in usual parametric amplification the system oscil- 
lates above threshold, that is for a gain g larger than f . 
This threshold is actually a consequence of the fact that 
the gain gT of parametric amplification becomes equal 
to the mechanical losses which are related to the damp- 
ing r. The damping T2 then vanishes and the noise of 
quadrature X2 becomes infinite. Above threshold the 
mirror oscillates at frequency f^f in phase with quadra- 
ture A2. The amplitude of oscillation depends both on 
the gain and on saturation mechanisms which take place 
in the parametric oscillation regime. In our case they are 
mainly related to the saturation of the intensity modula- 
tion of the auxiliary laser beam. 

The distribution of the mirror motion no longer ap- 
pears as a peak located at the center of phase-space 
(Xi = 0, X2 = 0). It rather looks like a squeezed dis- 
tribution centered at a non-zero position along the X2 
axis, the amount of squeezing depending on the distance 
to threshold. As the modulation of the feedback loop is 
at frequency 2Qm, the system is invariant by a 7r phase- 
shift and there are two possible positions for the center 
of the distribution, one for a positive value of (A^), the 
other one for a negative and opposite value — (A2). 

Figures Hni and lTTI show the experimental temporal evo- 
lutions and distributions obtained above threshold, re- 
spectively for a gain close to unity and for a larger gain. 
As expected the distribution in figure EH exhibits two 
peaks along the X2 axis, symmetrically located with re- 
spect to the center. The temporal evolution (top curve 
in figure ITU)! also shows a concentration of the mirror mo- 
tion around the two opposite positions. One can however 
note that the motion sometimes explores the neighbour- 
hood of the center of phase-space and can jump from one 
position to the other one. Delay between two jumps can 
be as small as a few seconds for a gain very close to unity 
and corresponds to a few minutes in the case of figure HIH 
Over the 10-minutes acquisition time of the distribution 
the system has actually made only one jump and it has 



stayed a longer time near the negative position than near 
the positive one. This explains why the negative peak is 
higher than the positive one. For a longer acquisition time 
one would obtain symmetric peaks. 

For a larger gain the system can stay a few hours with- 
out any jump. During the acquisition time the mirror 
then stays near one position and the distribution exhibits 
only one peak as shown in figure ITT1 




FIG. 10: Parametric oscillation of the mirror for a gain g > 1, 
in the {Xi,Xi) phase-space. Scale is 4 times larger than in 
figure |21 Top: temporal acquisition over 500 ms. Bottom: 
corresponding histogram acquired over 10 minutes. The dis- 
tribution appears as two opposite squeezed peaks with non- 
zero equal mean amplitudes. 
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Quadrature X! (10" 15 m) 




FIG. 11: Parametric oscillation of the mirror for a gain g much 
larger than 1, in the (Xi,X2) phase-space (same scale as in 
figure tTUt . Top: temporal acquisition over 500 ms. Bottom: 
corresponding histogram acquired over 10 minutes. 
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Parametric oscillation also squeezes the mirror motion 
along the X 1 quadrature. As usual the squeezing effect is 
larger close to threshold and the distribution goes back 
to a disk for large gain (figurc lTTl) . In contrast with para- 
metric oscillation in optics |32l l33| , the fluctuations are 
of the same order as the mean value of the oscillation, 
that is the distance between the two peaks in figure ITUI 
is of the same order as their widths. As a consequence 
one cannot use a linear approach to describe the evolu- 
tion of fluctuations and the distribution does not appear 
as an ellipse but as a more complex and distorted shape 
|34| . It is however possible to evaluate the variance AX^ 
and one obtains in the case of figure EI] a variance equal 
to 0.55 times the variance of the thermal noise at room 
temperature. 

Let us finally note that the amplitude of the peaks (dis- 
tance to the center of phase-space) depends on the gain. 
One gets a larger amplitude for a higher gain. There is 
however a limit related to the saturation of the feedback 
loop. For a very large gain the radiation pressure applied 
on the mirror is limited by the modulation capabilities 
of the laser beam intensity. One can easily estimate this 
limit. When the feedback loop saturates, the intensity 
modulation tends towards a square signal with a 100% 
modulation depth. The component of the radiation pres- 
sure force at frequency J7m is thus on the order of 

F rad (t)~ 2-- cos (fl M t), (35) 

C 7T 

where P is the light power. In this expression 2P/c repre- 
sents the constant force exerted by light without modula- 
tion and 4/7r is a correction factor due to the fact that the 
modulation corresponds to a square signal rather than to 
a sine one. The response of the mirror is a forced oscilla- 
tion with a non-zero mean amplitude (A 2 ) obtained by 
identifying the intrinsic damping force to the radiation 
pressure, 

MT%(X 2 )~2--. (36) 

C 7T 



For a 500 raW beam, this corresponds to a mean dis- 
placement of 6 10~ 15 m, in very good agreement with 
the one observed in figure ITT1 



VII. CONCLUSION 



We have presented an experiment of high-sensitivity 
displacement measurement, reconstructing the phase- 
space distribution of a mechanical oscillator with a 
thousand-fold increase in sensitivity upon what had pre- 
viously been reported. The sensitivity of our experimen- 
tal setup is currently limited to 1.7 10 -17 m by the damp- 
ing timescale of the oscillator, but it could easily be ex- 
tended to the attometer level by reducing the analysis 
bandwidth. 

Both feedback schemes studied in this paper are able 
to reduce the thermal noise of the mirror. Their effects 
however are different. In the case of cold damping the 
resulting state is still a thermal equilibrium but at a 
lower effective temperature, in principle down to a zero 
temperature for a very large gain |35j. One then gets a 
symmetric gaussian distribution in phase-space. In con- 
trast, parametric amplification only cools one quadra- 
ture of the mirror motion, the other one being heated. 
One then gets squeezed thermal states with an asymmet- 
ric gaussian distribution in phase-space. This situation 
is somewhat equivalent to the one which would be ob- 
tained with quadrature-dependent feedback loops where 
the zero-temperature quantum state of the mirror should 
be reached via squeezed thermal states [3(|. 

These results show that the observation of motion in 
phase-space gives a better understanding and control of 
optomechanical coupling, on the way to experimental 
demonstration of related quantum effects. 
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